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Аннотация 

Введение. Основным способом оценки несущей способности дорожных конструкций (представляющих собой 
слоистую среду, состоящую из материалов разной жесткости) является испытание установками ударного нагру- 
жения. Анализ результатов этих испытаний сводится к рассмотрению только пиковых значений вертикальных 
перемещений без учета изменения этих параметров в частотном и временном диапазоне. Вместе с тем именно 
изучение амплитудно-частотных характеристик перемещений может позволить решить ряд актуальных проблем, 
связанных, например, с приведением фактических значений вертикальных перемещений к значениям при про- 
ектных режимах работы конструкции. 

Целью исследования является анализ динамического отклика слоистых сред (дорожных конструкций) автомо- 
бильных дорог на ударное нагружение, направленный на приведение фактических параметров динамического 
отклика к проектным значениям, регламентируемым требованиями отечественных нормативных документов по 
проектированию дорог. 

Материалы и методы. Современные установки ударного нагружения позволяют воспроизводить нагрузку со 
временем импульса 0,03 с, в то время как проектные значения модулей упругости рассчитаны для 0,1 с. Для пре- 
одоления этого противоречия авторами было реализовано точное решение динамического уравнения Ламе для 
многослойного полупространства, и построен ряд амплитудно-частотных характеристик перемещений, на основе 
обработки которых была выведена упрощенная функциональная зависимость. 

Результаты исследования. Разработана методика построения амплитудно-частотных характеристик перемеще- 
ний на поверхности дорожной одежды при ударном воздействии, которая находит свое практические примене- 
ние. На основе точного решения динамического уравнения Ламе путем статистической обработки установлена 
упрощенная регрессионная зависимость для приведения вертикальных перемещений. Представлен пример рас- 
чета модулей упругости слоев дорожной конструкции с учетом их приведения в зависимости от времени прило- 
жения нагрузки. 

Обсуждение и заключения. В рамках исследования впервые получена методика приведения фактических значе- 
ний вертикальных перемещений к значениям при проектных режимах работы дорожной конструкции, целиком 
базирующаяся на расчетном анализе ее динамического НДС, в частности, амплитудно-частотной характеристике 
деформирования при нестационарном воздействии. Разработанный подход способен значительно повысить до- 
стоверность диагностики состояния автомобильных дорог, а также обеспечить единство измерений и результа- 
тов, получаемых при проектировании и диагностике дорог. 
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Введение. Одной из важнейших задач оценки прочности эксплуатируемых дорожных одежд на автомобиль- 
ных дорогах является корректный учет температуры и времени нагружения. Общеизвестно, что основным кри- 
терием, определяющим прочность дорожной одежды, является ее расчет на допускаемый упругий прогиб. На 
этапе проектирования модули упругости слоев асфальтобетона назначаются для расчетной температуры 10 °С и 
времени воздействия колеса расчетного автомобиля — 0,1 с. В условиях эксплуатации, когда осуществляется 
измерение прочности, фактическая температура и время воздействия от различных установок динамического 
нагружения могут существенно отличаться от этих нормативных значений, что несомненно вносит вклад в точ- 
ность и достоверность результатов измерений. 

Однако необходимо констатировать, что в последние годы исследований, касающихся учета влияния темпе- 
ратуры на свойства асфальтобетона, проведено крайне мало. Можно отметить исследования [1-3], позволившие 
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установить, что зависимость, связывающая общий модуль упругости дорожной одежды и температуру поверх- 
ности ее покрытия, имеет логарифмический характер и принимает следующий вид: 

Е =КхХ ШТ +В (1) 
где Т — температура связных слоев, °С; К, В — эмпирические коэффициенты. 

Коэффициенты К, В в свою очередь зависят от толщины связных слоев дорожной одежды. Помимо вышеука- 
занных исследований рекомендации по приведению модуля упругости асфальтобетонов к нормативным пара- 
метрам присутствуют в документе СТО АВТОДОР 10.1-2013 «Определение модулей упругости слоев эксплуа- 
тируемых дорожных конструкций с использованием установки ударного нагружения». В данном документе при- 
водятся рекомендуемые коэффициенты перехода от модуля упругости асфальтобетона, замеренного при факти- 
ческой температуре на поверхности покрытия и времени приложения нагрузки { = 0,03 с, к нормативной темпе- 
ратуре 10 °С и времени воздействия колеса расчетного автомобиля — 0,1 с. 

В зарубежной практике следует констатировать более широкий набор исследований, затрагивающих данную 
проблему. Так в исследованиях [4—6] указано, что модули упругости асфальтобетона, определенные в ходе натур- 
ных испытаний при фактической температуре, следует приводить к 20 °С либо к 25 °С. Механизм приведения в 
этих исследованиях также базируется на применении регрессионных зависимостей, связывающих модули упру- 
гости слоя с их температурой в диапазоне, как правило, от 0 до 40 °С. 

В работах [7—10] отмечается влияние времени приложения нагрузки на значение модуля упругости слоев ас- 
фальтобетона. Отмечается, что с увеличением времени приложения нагрузки снижается модуль упругости слоев 
асфальтобетона. В основе самого алгоритма приведения в представленных работах лежит соотнесение динами- 
ческого модуля упругости, зарегистрированного при частоте приложения нагрузки установки ударного нагруже- 
ния ЕУ/О, с мастер-кривыми динамического модуля упругости асфальтобетонных образцов, отобранных из по- 
крытия в точке ударного воздействия. В работе [11] приводятся сводные данные о принимаемых в различных 
источниках стандартных значениях частоты нагружения, воспроизводимой установками ударного нагружения 
Е\У/О (таблица 1). 

Таблица 1 
Данные о частоте нагружения установкой ударного нагружения ЕУ/О, принимаемые различными авторами 


Репрезентативная Частотно- 
Место исследований Источник частота нагрузки временное Тип Е\О 
РУ, Гц преобразование 

частота нагрузки определя- 
МоКОАР Супе, её а1., 2004 17,9 ется по временным данным Рупаез( 

тензометрического датчика 

Упепиа Тесь ен 
: . Гой|Н71, её а|., 2007 5,3 21 Оупаез( 
Тгапзромайоп шзййие 
Капзаз Зе Ошуегзйу Седа{а, её а1., 2010 25 ЕШ Рупаез( 
СЕРЕХ Маоз, её а1., 2012 15 +Ш КОАВ 

ны Аууайа, её а1., 2017 30 = Рупагез 


Аз$зос1аез, шс. (АКА) 


Как видно из представленной таблицы в исследованиях, выполненных на полигоне МпКоа4, а также в иссле- 
довательских организациях Капзаз 5‘ае Ошуегзиу, СЕОЕХ реализуется единый подход к определению частоты 
от приложения нагрузки установки ударного нагружения Е\/О, напрямую связывающий время импульса со вре- 
менем. Различие в значениях частоты можно связать с разным расчетным временем импульса воздействия уста- 
новки ударного нагружения. Иной подход применен в Утеииа ТесЬ Тгапзронаноп шзНвие. В соответствии с [12] 
авторами было установлено наличие сдвига на величину 1/2л по оси частот, что и отражено в приведенной рас- 
четной формуле. 

Таким образом, очевидна малая изученность вопроса выполнения временной и температурной корректировки 
результатов измерений модулей упругости в условиях эксплуатации в практике РФ и серьезная разрозненность 
применяемых подходов в международной практике. Это позволяет сделать вывод об актуальности изучения ме- 
ханизмов приведения значений модулей упругости слоев асфальтобетона по частоте и температуре к норматив- 
ным значениям. 

Материалы и методы. При проведении исследований была использована математическая модель динами- 
ческого напряженно-деформированного состояния многослойного полупространства. Данная модель базиру- 
ется на решениях, полученных в трудах [13—19]. Подробно постановка задачи и определяющие соотношения, 
реализуемые в математической модели, описаны в [20]. В рамках поставленной в данном исследовании задачи 
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анализируется динамическое напряженно-деформированное состояние многослойного полупространства под 
воздействием ударной нагрузки, распределенной по круговому отпечатку. Деформирование среды описыва- 
ется системой уравнений Ламе. Решение этих уравнений осуществляется с использованием преобразования 
Ханкеля в виде: 


Ме Ее 
094 (В, @) = | Л. ивУ (и, 5) .Х" лид 


(2) 
Г. п= 
где № — функция Бесселя; Ри") — ядро интегрального представления для многослойного полупространства; 
(п) 
Хх =: преобразования Ханкеля для напряжения на границах слоев (п =1, 2), определяемых из граничных усло- 


ВИЙ. 
Ядра интегрального преобразования для случая однородного полупространства имеют вид матриц-функций 
РИХи, 5): 


Р® (и, 5) = Р(и, <)- В" (и.0) 


(3) 
Виа —иехр(о' 5) — 01» ехр(о' <) 
’ оиехр(о“ 15) иехр(о'>5) 
(4) 
2 
В 0)=| $? 29а да в ПЕ. 
—2и0: -6> ‚Ав =с 2 4и 911012 52 =и тб (5) 


Преобразование Фурье от импульса хи) для однородного полупространства в форме преобразования 

Ханкеля имеет вид: 
Хи = ооо 
ь (6) 
где оо — давление на поверхность покрытия от воздействия установки ударного нагружения, Па; Р(&) — спектр 
нагружения. 
й (п) 

Вывод представлений Р9”, Х для многослойного полупространства осуществляется с применением ре- 
куррентных зависимостей, опирающихся на условие равенства перемещений на границах контакта слоев много- 
слойного полупространства. Общий вид зависимостей представлен в [21]. Применение интегрального преобра- 
зования Ханкеля позволяет осуществить построение амплитудно-частотной характеристики перемещений в лю- 
бой точке исследуемой среды. Амплитудно-частотная характеристика (АЧХ) при этом имеет вид, представлен- 
ный на рис. 1. 


200 400 600 800 1 000 
ми, рад/с 


Рис. 1. АЧХ перемещений на поверхности многослойного полупространства: 
9 7(\), мм — вертикальные перемещения; \, рад/с — круговая частота 
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Данный график построен для нагрузки величиной 700 кПа и времени приложения 0,03 с. Частотная область, 
в которой отслеживаются колебания поверхности покрытия, находится в диапазоне от 0 до 1500 рад/с (от 0 до 
238 Гц соответственно). Вертикальные перемещения изменяются в диапазоне от 0,43 мм при © = 0 рад/с (что 
соответствует статике) до 0,115 мм при вк = 1500 рад/с. Таким образом, в зависимости от частоты приложения 


нагрузки упругий прогиб изменяется в 3,74 раза. Для изучения влияния динамичности нагружения на величину 
вертикальных упругих деформаций был выполнен численный эксперимент для трех вариантов дорожных кон- 
струкций, отличающихся толщиной пакета асфальтобетонных слоев. Исследуемые конструктивы дорожных 
одежд и расчетные параметры представлены в таблице 2. 


Таблица 2 
Входные данные для моделирования 
Конструкция дорожной одежды Толщина слоев, см д и Е 
коэффициент Пуассона демпфирования 
Асфальтобетон 10-30 2000/0,35 0,08 
Основание, не обработанное вяжущим 40 200/0,35 
Грунт земляного полотна — 41/0,35 


Модуль упругости слоя асфальтобетона принимался в виде средневзвешенного значения, коэффициент 
демпфирования слоя асфальтобетона задавался в соответствии с [22]. Результаты моделирования представлены 


в таблице 3. 
Таблица 3 
Результаты моделирования амплитудно-частотной характеристики деформирования на поверхности покрытия 
Конструкция 1 Конструкция 2 Конструкция 3 
Е й 1 и Кпер Е р 1 и Кпер Е й 1 и Кпер 

0 | 0,562 | 1,070476 0 | 0,59 | 1432039 0 | 0,48 | 1.142857 
30 | 0,5 | 1.106195 30 | 0,462 | 1,121359 30 | 0,385 | 1.166667 
0 60 | 0,395 | 1,09116 0 60 | 0,362 | 1.113846 0 60 | 0,32 | 1,185185 
120 | 0,255 | 1,20283 120 | 0,248 | 1,1698 120 | 0,235 | 1,236842 
180 | 0,175 | 1,129032 180 | 0,175 | 1.166667 180 | 0,17 | 1.214286 

0 | 0,525 1 0 | 0,53 | 1.286408 0 | 0,42 1 

30 | 0,452 1 30 | 0,412 1 30 | 0,33 1 

10 60 | 0,362 1 10 60 | 0,325 1 10 60 | 0,27 1 

120 | 0,212 1 120 || 0.212 1 120 | 0,19 1 

180 | 0,155 1 180 [0,15 1 180 [0,14 1 
0 | 0,42 0,8 0 | 0,412 1 0 | 0:3 | 0,714286 
30 | 0,355 | 0,785398 30 | 0,305 | 0,740291 30 | 0,23 | 0,69697 
25 60 | 0.251 | 0,69337 25 60 | 0.23 | 0.707692 25 60 | 0.18 | 0.666667 
120 | 0,151 | 0,712264 120 | 0,12 | 0,566038 120 | 0,12 | 0,631579 
180 [0,11 | 0.709677 180 | 0,095 | 0,633333 180 | 0,08 | 0,571429 
0 | 0:3 | 0,571429 0 | 0,305 | 0,740291 0 | 0,215 | 0,511905 
10 | 30 | 0,25 | 0,553097 20 | 30 | 0,205 | 0.497573 30 |_ 30 | 0.145 | 0,439394 
50 60 | 0,175 | 0,483425 50 60 | 0,145 | 0,446154 50 60 | 0,11 | 0,407407 
120 | 0,095 | 0,448113 120 | 0,08 | 0,377358 120 | 0,065 | 0,342105 
180 | 0,062 0,4 180 | 0,05 | 0,333333 180 | 0,04 | 0,285714 
0 |0 0,4 0 | 0,212 | 0,514563 0 | 0,14 | 0,333333 
30 | 0,162 | 0,358407 30 | 0,0125 | 0,03034 30 | 0,085 | 0,257576 
100 60 | 0,1 | 0,276243 100 60 | 0,075 | 0,230769 100 60 | 0,06 | 0,222222 
120 | 0,05 | 0,235849 120 | 0,04 | 0,188679 120 | 0,03 | 0,157895 
180 | 0,037 | 0,23871 180 | 0,025 | 0,166667 180 | 0,02 | 0,142857 
0 | 0,162 | 0,308571 0 | 0,152 | 0,368932 0 | 0/11 | 0,261905 
30 | 0,11 | 0,243363 30 | 0,075 | 0,182039 30 | 0,05 | 0,151515 
180 60 | 0,06 | 0,165746 180 60 | 0,048 | 0,147692 180 60 | 0,03 | 0.111111 
120 | 0,047 | 0,221698 120 | 0,25 | 1/179245 120 | 0,02 | 0,105263 

180 | 0 0 180 0 0 180 | 0 0 
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Моделирование осуществлялось для набора частот 0, 10, 25, 50, 100, 180 Гц и для точек наблюдения, находя- 
щихся на расстоянии от точки приложения нагрузки (1) — 0, 30, 60, 120, 180 см. В качестве стандартной частоты 
принималась частота 10 Гц, и для нее же считался коэффициент перехода Кер. По итогам статистической обра- 
ботки полученных результатов была построена корреляционная зависимость К (р [, й) между коэффициентом пе- 
рехода от фактической частоты нагружения к расчетной частоте 10 Гц, имеющая следующий вид: 


К(Е,1, В) = 1,305738 — 0,014451} + 0,000052{2 — 0,066578 т(й) — 0,0004671 (7) 


где Г — частота колебаний, Гц; й — толщина пакета слоев асфальтобетона; [ — расстояние от точки ударного 
нагружения. 

Статистическая обработка осуществлялась в программном комплексе З‘айзйса. Коэффициент детерминации 
данной зависимости составил К2 = 0,97. Расчет и приведение вертикального перемещения к нормативной частоте 
и времени нагружения осуществляется по формуле: 


Онорм В И 25гц х к(р, 1, в) (8) 
Результаты исследования. В качестве примера реализации полученной зависимости был осуществлен рас- 
чет и приведение к расчетной частоте чаши прогибов и модулей упругости слоев дорожной одежды на эксплуа- 
тируемом участке автомобильной дороги М1 «Беларусь». 
Конструкция дорожной одежды, для которой осуществлялся расчет, представлена в таблице 4. 


Таблица 4 
Конструкция дорожной одежды, принимаемая при моделировании 
р Проектное значение модуля упругости, 
Слой Толщина, см 
МПа 

Асфальтобетон 19 2500 
Основание 33 200 
Дополнительный слой основания 50 120 
Грунт земляного полотна 41 


Фактическая чаша прогибов, зарегистрированная при фактическом времени контактного воздействия 0,04 с 
(что соответствует 25 Гц), и чаша прогибов, откорректированная в соответствии с зависимостью (7), представ- 
лены в таблице 5. 

Таблица 5 


Результаты моделирования амплитудно-частотной характеристики деформирования 
на поверхности покрытия 


п Датчик | Датчик 1 | Датчик | | Датчик! | Датчик! | Датчик | | Датчик | | Датчик | | Датчик 1 
арамет 
Г ыы 1 (0) (200) (300) (600) (900) (1200) (1500) (1800) (2100) 
Прогиб, мм, 
при частоте 0,380 0,326 0,291 0,245 0,210 0,145 0,102 0,072 0,053 
25 Гц 
Коэффициент 
ПО бар. |. м 0,767 0,753 0,739 0,725 0,711 0,697 0,683 
частоте 
к, Ьв) 
Прогиб, мм, 
при частоте 0,487 0,423 0,379 0,325 0,284 0,200 0,143 0,103 0,078 
10 Гц 
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В,м 


— Расчет 
— Экспер 


Рис. 2. Расчетная и экспериментальная чаша прогибов до корректировки 


Рис. 3. Расчетная и экспериментальная чаша прогибов после корректировки 


Таблица 6 


Результаты определения модулей упругости конструктивных слоев дорожных одежд от фактического 


времени приложения нагрузки 


Откорректированные модули 


Откорректированные модули 


Грунт земляного полотна 


ыы (= 0,03 с) (@(=0,1 с) 
Асфальтобетон 3200 3014 
Основание 348 190 
Дополнительный слой основания 62 55 
170 95 


Как ВИДНО ИЗ таблицы 2, приведение модулей упругости слоев дорожной одежды по времени в соответствии 
с полученными зависимостями затрагивает все материалы конструктивных слоев, так как для многослойной си- 
стемы изменение жесткости отдельных слоев влечет за собой изменения в форме АЧХ перемещений, что про- 
слеживается в данных, представленных в таблице 1. Полученная зависимость позволяет осуществить расчет 
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общего модуля упругости на поверхности дорожной конструкции для времени приложения нагрузки / = 0,03 с, 
и {= 0,1 с. В этом случае общий модуль упругости рассчитывается как: 
ро(а-и?) 

Еобщ —= о (9) 
где Еоыц — общий модуль упругости на поверхности дорожной конструкции, МПа; р — давление, кПа; р — диа- 
метр отпечатка, м; и — коэффициент Пуассона; 4о — упругий прогиб в центре приложения нагрузки, мм. 

В этом случае значение общего модуля упругости при времени воздействия {1 = 0,03 с составит 575 МПа, а 
значение общего модуля упругости при времени воздействия {1 = 0,1 с — 449 МПа. Подобное приведение необ- 
ходимо, когда речь идет о сопоставлении фактического общего модуля упругости с проектным значением, опре- 
деленным по нормативным документам. 

Обсуждение и заключения. Разработана методика построения амплитудно-частотных характеристик пере- 
мещений на поверхности дорожной конструкции при ударном воздействии. 

Выполнено численное моделирование изменения амплитудно-частотной характеристики перемещений на по- 
верхности дорожной конструкции для различных сочетаний толщины асфальтобетонных слоев, расстояния от 
точки ударного воздействия и частоты колебаний. 

На основе точного решения динамического уравнения Ламе путем статической обработки установлена упро- 
щенная регрессионная зависимость для приведения вертикальных перемещений, зарегистрированных при ча- 
стоте 25 Гц (соответствующей частоте нагружения установки Е\\/О), к нормативному времени нагружения 0.1 с 
(соответствующему нормативным документам). Коэффициент детерминации составил Е2 = 0,97. 

Представлен пример расчета модулей упругости слоев дорожной конструкции с учетом их приведения в за- 
висимости от времени приложения нагрузки. Установлено, что при приведении модулей упругости слоев дорож- 
ной одежды от 25 Гц к 10 Гц происходит их снижение, проявляющееся для всех материалов слоев. 
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